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Quantum anomalies are the breaking of a classical symmetry by quantum fluctuations. They 
dictate how physical systems of diverse nature, ranging from fundamental particles to crystalline 
materials, respond topologically to external perturbations, insensitive to local details. The anomaly 
paradigm was triggered by the discovery of the chiral anomaly that contributes to the decay of pions 
into photons and influences the motion of superfluid vortices in 3 He-A. In the solid state, it also 
fundamentally affects the properties of topological Weyl and Dirac semimetals, recently realized 
experimentally in TaAs, NaaBi, Cd3As2, and ZrTes- In this work we propose that the most identi¬ 
fying consequence of the chiral anomaly, the charge density imbalance between fermions of different 
chirality induced by non-orthogonal electric and magnetic fields, can be directly observed in these 
materials with the existing technology of photoemission spectroscopy. With angle resolution, the 
chiral anomaly is identified by a characteristic note-shaped pattern of the emission spectra, origi¬ 
nating from the imbalanced occupation of the bulk states and a previously unreported momentum 
dependent energy shift of the surface state Fermi arcs. We further demonstrate that the chiral 
anomaly likewise leaves an imprint in angle averaged emission spectra, facilitating its experimental 
detection. Thereby, our work provides essential theoretical input to foster the direct visualization 
of the chiral anomaly in condensed matter. 

Quantum anomalies occur when a classically preserved 
symmetry is broken upon quantization 1 ! The chiral 
anomaly, the first known example, refers to the non¬ 
conservation of a chiral current — a current imbalance 
between two distinct species of chiral fermions. Its dis¬ 
covery resolved the discrepancy between the measured 
and calculated decay rate of the neutral pion into two 
photon^, and in condensed matter it was first identi¬ 
fied as being responsible for a transverse force acting on 
vortices in superfluid 3 He-A®. In the solid state, Weyl 
and Dirac semimetals—TaAs was recently theoretically 
predicted ^ 5 an j experimentally identified using angle- 
resolved photoemission spectroscopy^ (ARPES) as be¬ 
ing a Weyl semimetal, with Na3BP®l, Cd3 AsJEHIHI and 
ZrTes 19 likewise identified as Dirac semimetals—are de¬ 
scribed at low energies by chiral fermions that realize a 
chiral ano maly and host topological surface states with 
Fermi arcs 20 21 . In this work we theoretically explore the 
spectroscopic signatures of the chiral anomaly in these 
materials, with special emphasis on its effect on the Fermi 
arcs. We reveal an identifying note-shaped structure ob¬ 
tained in ARPES and a characteristic peaked structure 
of angular averaged photoemission spectra that allow for 
an experimental visualization of the chiral anomaly. 

Weyl semimetals are a topological state of matter in 
which the conduction and valence bands touch and lin¬ 
early disperse around pairs of Weyl nodes 20 21 . Each 
node has a definite left or right handed chirality pro¬ 
viding a quantum number analogous to the valley de¬ 
gree of freedom in graphene 22 . Dirac semimetals can be 
thought of as two superimposed copies of Weyl semimet¬ 
als with the degeneracy pr otect ed by a crystal symme¬ 
try from opening up a gapP 5 23 25 . Similar to topologi¬ 
cal insulators and their metallic surface, Dirac and Weyl 
semimetals host protected surface stated^ that uncon¬ 


ventionally only exist for a restricted range of crystal 
momenta, thereby forming a Fermi arc connecting a pair 
of Weyl points with opposite chirality 2 ®^. Fermi arcs 
have been identified with ARPES in both TaAs ^ 7 and 
Na 3 BP. 

The chiral fermions describing the low energy degrees 
of freedom of Dirac and Weyl semimetals exhibit the chi¬ 
ral anomal}®®: while the sum of left and right handed 
fermions is necessarily conserved, their difference, the 
chiral density, does not have to be, even if classically 
it should. In fact, non-orthogonal magnetic and elec¬ 
tric fields pump left handed fermions into right handed, 
or vice verssPP Disorder induced internode scatter¬ 
ing eventually counterbalances the pumping, leading to 
a non-equilibrium steady state with a nonzero chiral 
density. The chiral anomaly is theor etical ly predicted 
to result in negative magnetoresist anc^^ ^ (experimen¬ 
tally observed in ZrTes 19 and TaA d 37 * 3 ^ ) , loca l 39 * 40 -! and 
no n-locaP - 1 transport phenomena, chiral optical activ- 
itp 2 42 * 43 * and rotation induced coolinJHI. 

Photoemission spectroscopy (PES) has previously 
been overlooked as a diagnostic tool of the chiral anomaly 
despite its importance in observing topological semimet¬ 
als. A possible reason for this is that the finite magnetic 
field required to observe the chiral anomaly may compli¬ 
cate the disentangling of electron trajectories needed for 
angular resolution. Setting aside this obstacle (we return 
to it at the end of this work where we argue that it can 
be overcome) we identify the main spectral signatures of 
the chiral anomaly and its observable effect on the Fermi 
arcs. First, the bulk spectrum is determined by the dif¬ 
ferently occupied Weyl nodes distinguished by their chi¬ 
ralities; second, the bulk chiral imbalance tilts the Fermi 
arcs, which then appear at fixed energy as finite segments 
stemming from the bulk Fermi surface. Together, these 
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FIG. 1. Visualization of the chiral anomaly in Dirac and Weyl semimetals, (a) Low-energy spectrum of a Weyl 
semimetal film with two bulk Weyl nodes of different chirality separated in momentum space. The grey plane represents the 
surface state at the film’s top surface, occupied up to the equilibrium chemical potential /r eq . Applying external magnetic and 
electric fields that satisfy E'B / 0 results in a steady state with left and right cone chemical potentials \il ^ f^R, linearly 
interpolated by a tilted Fermi arc. (b) Two constant energy cuts (A and B) through the band structure, with occupied and 
empty surface states depicted by solid light blue and white dashed lines respectively. The occupation at these cuts shows a 
characteristic blue note-shaped pattern depicted in the lower panel, (c) Dirac semimetals host pairs of Weyl cones, each pair 
with fixed isospin (t or and both left and right chiralities, that respond to the chiral anomaly in the opposite way. Two edge 
states with opposite velocities (light red and light blue planes), appear at each boundary of the Dirac semimetal. Scattering 
processes within and between cones with scattering times r c , t v and Ti are depicted by arrows, (d) The two pairs of Weyl 
nodes in (c) together comprise a pair of Dirac nodes. At a fixed energy cuts (C and D) between jiL and /ir , both bulk nodes 
are occupied while the surface states are only partially occupied. The total occupation in these planes describes two facing 
note-shaped patterns, illustrated in the bottom panel. 


two features form a distinguishing note-shaped photoe¬ 
mission pattern that we argue are within reach of current 
experimental state of the art. We further calculate angle 
integrated PES, which does not suffer from the magnetic 
field complications of its ARPES relative, and show that 
distinct signatures of the chiral anomaly survive. Overall, 
our results supply essential theoretical input that render 
photoemission spectroscopy a viable probe to visualize 
the chiral anomaly in Dirac and Weyl semimetals. 

To support our conclusions, we start by discussing the 
main features of the band structure of Weyl semimet¬ 
als, shown in Fig. EE Their low-energy bulk spectrum 
consists of an even number of band touching points of 
left (L) and right (R) handed chirality that are sepa¬ 
rated in energy and momentum by breaking inversion or 
time-reversal symmetry. Close to the Weyl points the 
energy dispersion is approximately linear and described 
by a Weyl Hamiltonian 

Hl/r = ±hva- ■ k, (1) 

with k the momentum, a the vector of Pauli matrices 
and v the Fermi velocity. Its eigenstates have a spin that 
points either radially away from or into the Weyl point, 
depending on its chirality. Each Weyl node therefore acts 
as a monopole of Berry flux in momentum spacd^ 45 . In 


analogy with magnetic monopoles, a Dirac string neces¬ 
sarily emanates from the monopoles connecting a pair 
with opposite chirality—. Any two dimensional plane in 
momentum space, spanned say by the two momenta k x 
and k y at a fixed k z , that crosses the Dirac string an 
odd number of times and does not contain a Weyl point 
defines a topologically nontrivial gapped band structure. 
Therefore, if the system is made into a film that is finite 
in a real space direction conjugate to a momentum in that 
plane, for example either the x- or ^/-direction, a chiral 
surface state is obtained at the corresponding surfaces 26 . 
Since k z remains a good quantum number we can repeat 
this arguments for different planes at different k z . Only 
those crossing the Dirac string odd times have a surface 
state, which therefore only exists for certain values of the 
k z momenta. The surface state dispersion is depicted as 
the grey shaded plane leaning on the two Weyl nodes 
in Fig. Ef (the opposite surface provides an analogous 
surface plane with opposite velocity that is not shown). 
The separator between occupied and unoccupied surface 
states is an arc — the Fermi arc. 

In the presence of non-orthogonal external electric (E) 
and magnetic (B) fields, the chiral anomaly leads to a 
non-conservation of the left and right handed electron 
densities. This is expressed by the two coupled continuity 
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FIG. 2. ARPES signatures of the chiral anomaly in Weyl semimetals. Numerically computed ARPES spectra for a 
doped Weyl semimetal film with L = 2000 layers in (a) equilibrium and (b) and (c) the chiral anomaly induced steady state. 
The parameters are such that the equilibrium chemical potential is /x eq = 0.02 v while in non-equilibrium the left and right 
cones are filled up to rl = 0.025 v and rr — 0.008 v respectively. The lower panels in (a) and (b) show the momentum resolved 
spectra at a fixed energy, located as schematically shown with grey lines in the upper panels. In equilibrium (a), bulk cones 
and the surface state are observed, while the chiral anomaly (b) results in the disappearance of one cone and the emerging of 
the characteristic note-shaped pattern between rl and rr. The panel (c) shows the ARPES spectrum in non-equilibrium at 
a fixed k y — 0, depicted in the top panel. The upper plot in the lower panel shows the total ARPES spectra, demonstrating 
the cone-like structure from the bulk and the flat surface state; the lower plot displays the intensity difference A I = I ca — / eq 
between the out-of-equilibrium and the equilibrium states. The dashed lines mark the chemical potentials Rl, Rr and R eq . 
These plots were obtained for v = v z , t = 0.5 v, e = 6t and b z = 0.9 v with the temperature set to T — O.OOlu, which 
corresponds to T « 5 K for a typical Fermi velocity of v = 0.45 eV (see methods and supplementary material for details). 


equations 

3 

ed t n L/R + V-j L/R = ±j-(n R -n L ), (2) 

where n L / R is the density of left and right handed 
fermions measured from the Weyl point and Jl/r their 
current density. The first term on the right hand side is 
the anomaly contribution, which has a different sign for 
the two chiralities; the second term represents the inter¬ 
valley scattering with rate r~ l . At long times, a steady 
state with occupation difference between the two chiral¬ 
ities is obtained. The continuity equations © and par¬ 
ticle conservation then define the left and right handed 
chemical potentials 


I^l/r — 


MegT 2^ 3 e 2 r„ E B 


1/3 


(3) 


where R eq is the equilibrium chemical potential and we 
have used that n R / R = n \/ R /(6 7r 2 h 3 v 3 ) for three dimen¬ 
sional Weyl fermions. In defining Rl/r we assume that 
the equilibration 47 within a node with intra-valley relax¬ 
ation time r c is much shorter than the inter-valley equili¬ 
bration time r v 41 [ This similarly suggests that relaxation 
along the surfaces is dominated by small momentum scat¬ 
tering such that the Fermi arc linearly interpolates in 
momentum space between r r and r r , analogous to the 


voltage drop along an ohmic wire, leading to the steady- 
state occupation shown in the bottom of Fig. [l^i. At a 
fixed energy between r r and rl this results in the char¬ 
acteristic note-shaped occupation schematically shown in 
Fig[l]3 for the two constant energy planes denoted cut A 
and cut B. 

Dirac semimetals can be understood as two copies 
of Weyl semimetals, with each Dirac node composed 
of two Weyl fermions of opposite chirality. The two 
copies are distinguished by their different total angular 
momentum!^, which can be captured by an isospin quan¬ 
tum number that we denote with t and As long as the 
crystal symmetry is not broken, the isospin remains a 
good quantum number. This has two important conse¬ 
quences: one, the two chiral fermions comprising a Dirac 
node are decoupled and therefore a gap does not open; 
two, the Weyl nodes still act as monopoles in momentum 
space and their Dirac string connects monopoles with op¬ 
posite chirality but the same isospin. Dirac semimetals 
therefore have two Fermi arcs with opposite velocity on 
each surface. 

The two pairs of Weyl fermions in the Dirac semimetal 
are oppositely affected by the chiral anomaly. External 
and non-orthogonal E and B fields shift the occupation 
in one Dirac cone to higher energies for one isospin and 
to lower for the other, and oppositely in the other Dirac 
cone, as schematically shown in Fig. (TJ:. For this steady 























































4 



FIG. 3. ARPES signatures of the chiral anomaly in Dirac semimetals. Numerically computed ARPES spectra for 
a doped Dirac semimetal film of L = 2000 layers in (a) equilibrium and (b) and (c) the presence of E • B ^ 0. The spectra 
in (a) and (b) are the momentum resolved spectra at the fixed energies schematically shown at the top as gray lines. The 
equilibrium chemical potential is set to /z e q = 0.020 v while the left and right out of equilibrium chemical potentials are chosen 
to be )il = 0.025 v and /jlr = 0.008 v. While there is a slight bulk Dirac cone intensity reduction from the equilibrium to the 
non-equilibrium situation, a stark qualitative difference is observed in the surface states that leads to the characteristic double 
note-shaped pattern. In (c) we plot the intensity difference A I — / ca — / eq for a fixed momentum k z = ±0.5, which reveals the 
linear dispersion of the edge state. These plots were obtained for Mo = — 0.2 v, Mi = —0.25u, M 2 = —0.75u and a temperature 
of T — 0.001 v. With v — 0.45 eV for NasBi, the temperature corresponds to T « 5K and the induced chemical potential 
difference to 8meV. 


state to be realized, intra-valley relaxation at fixed chi¬ 
rality must be larger than both the inter-valley relax¬ 
ation, r” 1 r” 1 , and the intra-valley relaxation be¬ 

tween isospins, r c _1 >> r- _1 (the different relaxation pro¬ 
cesses are depicted with arrows in Fig. []J). Both condi¬ 
tions are estimated to be satisfied in Dirac seniimetals 41 . 
Moreover, high mobility materials such as graphene sat¬ 
isfy the condition r c _1 r” 1 , which allows for the exper¬ 
imental observation of a chemical potential imbalance 
and suggests that it is also met in CdsAs 2 . 

The total occupation of a given Dirac cone is a super¬ 
position of both isospins, and therefore the bulk occupa¬ 
tion in the chiral anomaly induced steady state is qual¬ 
itatively the same as in its absence: two circular disks. 
In contrast, the tilt of the Fermi arcs and the resulting 
partial occupation at a fixed energy, see Fig. [l]i, leads to 
a qualitatively new signature in the form of two facing 
note-shaped patterns (familiar from the Weyl semimetal 
case). This key property, which allows photoemission of 
Dirac semimetals to show evidence of the chiral anomaly, 
is a central result of our work. 

We establish by a numerical computation that these 
note-shaped patterns are indeed directly manifested in 
photoemission spectra. Our simulated data is obtained 
with exact diagonalization of tight binding models of 
Weyl and Dirac semimetals, as explained in the meth¬ 
ods section. In Fig. [2^i we plot the momentum resolved 
ARPES spectra at various fixed energies for a doped Weyl 
semimetal in equilibrium. Two bulk cones and one sur¬ 


face state at the probed surface are clearly seen; the sur¬ 
face state localized at the opposite surface is not visible 
due to the finite penetration depth of the incoming pho¬ 
ton. In the presence of external fields E-B^0, the note¬ 
shaped pattern of the occupied states is seen in Fig. 
for energies between and [ 1 r. An alternative way to 
illustrate the steady state occupation with /ir 7 ^ /jll is 
through a cut in momentum space at a fixed k y = 0 as 
provided in Fig. H- The pumping of charge between the 
different chiralities is clearly revealed in the intensity dif¬ 
ference between the equilibrium J eq and chiral anomaly 
induced steady state J ca , shown in the lower panel of 
Fig.[2j;. 

The numerically computed ARPES spectra for a Dirac 
semimetal, shown in Fig. [3} s imilarly reveal the essential 
features discussed in Fig7|TJ In equilibrium two bulk 
cones and two counter-propagating edge states at the 
same surface are seen in Fig. [3^t, similar to the exper¬ 
imental observations for NasBP 1 . The non-equilibrium 
spectra in the presence of electric and magnetic fields are 
qualitatively different; two copies of the note-shaped pat¬ 
tern clearly reveal the chiral anomaly in Dirac semimet¬ 
als. To further highlight this feature we plot, as for 
the Weyl semimetal, the intensity difference between 
the equilibrium and non-equilibrium in Fig. [ 3 J 3 for fixed 
k z = ±0.5. The qualitatively most notable feature is the 
partial occupation of the surface state that results in the 
stem of the note. The bulk, though, is not entirely insen¬ 
sitive to the chiral anomaly: first, states are occupied up 
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to an energy max(/i£, rr) that is higher than in the equi¬ 
librium situation; second, only one isospin band is filled 
between rr and rr resulting in a decreased intensity. 

From the chiral anomaly equation © we estimate the 
induced chiral chemical potential difference Sr = rl~ HR 
to be within experimental state of the art. For the dop¬ 
ing levels of NasBi, an electric field strength of 10 4 V m -1 
and a magnetic field of 1 mT gives Sr of the order of 
10 meV, which is well within ARPES resolution. Re¬ 
markably, magnetic fields as small as 6 rT can still 
achieve observable Sr ~ 2.7 meV for NasBi and can reach 
Sr ~ 5.4 meV in CdsAs 2 (details of the experimental val¬ 
ues used to compute these estimates for different mate¬ 
rials as well as ARPES spectra specific to NasBi with 
curved Fermi arcs, are given in the supplementary ma¬ 
terial). The fact that experimental observation of Fermi 
arcs in equilibrium have already been reported also bodes 
wel^M 

To achieve momentum resolution in ARPES it is nec¬ 
essary to correlate the angle at which an electron is de¬ 
tected to its initial momentum. In free space, electrons 
move in straight lines making this task straightforward, 
and a constant electric field does not overly complicate it. 
A magnetic field turns the electron trajectories into spi¬ 
rals and may, depending on the distance to the detector 
and size of the magnetic field, make momentum resolu¬ 
tion difficult. However, since the magnetic field strength 
needed for an observable effect is rather small, correcting 
induced deviations is plausible 49 . Alternatively, a larger 
uniform magnetic field finite within the sample could in 
principle be engineered with a ferromagnetic material in 
the picture frame geometry: a closed magnetic circuit 
that minimizes stray fields outside the sample that has 
been used to experimentally study ferromagnetic metals 
via ARPE^ni 

(see supplementary material for details). 
Similarly, if either field can be turned off fast enough 
(faster than the inter-valley relaxation) a pump-probe 
setup could observe the non-equilibrium steady state and 
its equilibration, which would allow for a direct measure¬ 
ment of the inter-valley relaxation time in addition to 
visualizing the chiral anomaly. While none of these are 
simple tasks, we believe that the rewards are significant 
enough that the experimental challenge will be met. 

In the absence of momentum resolution the angular av¬ 
eraged but energy resolved photoemission spectroscopy 
likewise contains direct signatures of the chiral anomaly. 
In Fig. [4] we plot the numerically computed PES spectra 
as a function of energy for both Weyl and Dirac semimet¬ 
als. In equilibrium it has a single step that is smeared 
by temperature. In the non-equilibrium steady state oc¬ 
cupation is shifted from lower energies to higher, such 
that the bulk spectra would have a double step profile, 
with one step dX rl and the other at Rr. In the total 
spectra the two steps are hard to see since the surface 
states contribute significantly to smoothen the profile. 
The intensity difference A I = J eq — J ca instead shows 
a characteristic peak-dip structure for low temperatures, 
which reflects the chiral anomaly pumping of fermions of 



FIG. 4. Chiral anomaly in PES spectra for Dirac and 
Weyl semimetals. Numerically computed PES intensities 
for Weyl and Dirac semimetals for two different temperatures. 
In equilibrium the intensity / eq presents a single step around 
fi eq , that is smeared out by increasing temperature and re¬ 
sembles the occupation of the cones. In non-equilibrium the 
intensity / ca shows a double-step-profile with steps at rr and 
rl- Both Rr/l include the temperature dependence of the 
chiral anomal}! 2 19 . The difference A I — / eq — / ca highlights 
the effect of the chiral anomaly by a characteristic peak-dip 
structure for low temperatures, resulting from the shift in oc¬ 
cupation of the two cones, that evolves into a single peak 
as temperature is increased. The parameters used to obtain 
these plots are /x eq = 0.020 v, rl — 0.025 v, rr — 0.008 v at 
4 K and rl — 0.024 v, rl — 0.005 v at 15 K. The temperatures 
are chosen to be T = 8 • 10 -4 v and T = 2.9 • 10 -3 v, corre¬ 
sponding to 4K and 15 K for a typical value of v = 0.45 eV. 
The rest of parameters are chosen as in Fig. [2] and [3] for the 
upper and lower panel respectively. 


one chirality into the other, evolving into a single peak 
as temperature is increased. 

In conclusion, we have promoted photoemission as an 
ideal tool to experimentally visualize the chiral anomaly 
in the solid state. The main effect that allows this 
is the tilt of the Fermi arc in the presence of non- 
orthogonal electric and magnetic fields that, through the 
chiral anomaly, pump electrons from one chirality to the 
other. As a consequence, the surface state occupation at 
a fixed energy no longer connects bulk nodes but rather 
terminates in between, resulting in a qualitatively dis¬ 
tinct note-shaped photoemission pattern identifying the 
chiral anomaly. We argued, by estimating from experi¬ 
mentally available parameters the relevant chiral chemi¬ 
cal potential difference to be about 10 meV, that a direct 
visualization of the chiral anomaly is within the state 
of the art. Were it to be realized experimentally in the 
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way proposed here, a revealing light would be shed on ky- and isospin-dependent. For the surface states, /jl^ 
the relation between two fundamental concepts: quan- depends linearly on k z . 
turn anomalies and topological states of matter. 


METHODS 
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^wsm = Hq + Hi, (4a) 

Hq — 2 u(sin ky cr x sin k x Gy) t~x H - 2 v z sin k z Ty H - A^\tT X: 

(4b) 

Hi = b 0 a z t v + b • (~<j x t x , a y r x , a z ). (4c) 


Appendix A: Aspects of experimental 
implementation 

1. Relaxation rates 


The first term models a three-dimensional topological in¬ 
sulator 51 . The Pauli matrices a and r respectively act 
in spin and particle-hole space. The velocity v in the x- 
and //-directions differs from that in the ^-direction (v z ) 
consistent with experimentally relevant materials. The 
mass Mk = e — 2 1 JT cos At e = 6t a degenerate 
three-dimensional Dirac cone is obtained at k = 0. The 
second terrrPS breaks inversion with a nonzero bo and 
time-reversal symmetry with a finite b. In our simula¬ 
tion we take bo = 0 and fix the direction of b = (0, 0, b z ) 
such that the Dirac cone is split into two Weyl cones 
at k z ~ The low-energy spectrum of both Dirac 

semimetals NasBi^ and CdaAsJ^ is modeled by the 
Hamiltonian 

Mdsm = M k a z + 2 v (sin k x cr x r z - sin k y <r y ), (5) 

which can be understood as two copies of two-band Weyl 
semimetal Hamiltonians. The mass term Mk = Mo + 

2 Mi (cos k z — l)+2 M 2 (cos k x +cos k y — 2) sets the velocity 
in k z direction to v z = 2 \/Mo Mi and the Dirac cones to 
be at k c z = ±a/M 0 /Mi. In this simple model the surface 
states are dispersionless in the /^-direction. 

To model a film we take the system finite in the direc¬ 
tion orthogonal to the separation of the cones in momen¬ 
tum space, which is the surface where the edge states are 
maximally visible. The ARPES spectrum is given by an 
integral over the local density of states 

7(k||,w)= / da: ^ \^nM\\ ( x )| 2 ^( e n, k|l ), 

** n 

(6) 

where x is the coordinate in the finite direction of the 
film, w(x) = exp(—x/£) is a weight function modelling 
the incoming light’s intensity decay with depth into the 
sample with decay length £, and / is the Fermi-Dirac 
distribution. The wavefunctions VVky are the eigenfunc¬ 
tions of the Hamiltonians and e n5 k|| the corresponding 
eigenvalues, which depend on the momentum ky parallel 
to the surface. To model the chiral anomaly induced 
steady state, the chemical potential fi k y is taken to be 


The experimental feasibility of detecting the chiral 
anomaly with photoemission spectroscopy relies firstly on 
the correct hierarchy of the different relaxation rates in¬ 
volved. As discussed in the main text, for Weyl semimet¬ 
als the intra-valley relaxation rate at fixed chirality must 
be faster than the inter-valley relaxation rate, r” 1 
r” 1 . For Dirac semimetals it is also required that the 
relaxation rate between the two isospins forming each 
Dirac node must be shorter than the intra-valley relax¬ 
ation, Ti r c . The intra-valley relaxation relaxation 
rate r c can be deduced from the experimental values for 
the carrier mobilities fi e given in the Table [I] using that 


a 

Me = — 


en 


e vf 
h kp Tc ’ 


(SI) 


where a denotes the DC conductivity, e is the electric 
charge, and n is the carrier density. A theoretical esti¬ 
mate of the ratio r v /r c can be calculated in the first Born 
approximation following Ref. EH rendering the values in¬ 
cluded in the Table [T] that justify the assumptions used 
in the main text. 

Experimentally, r v can be determined via non-local 
transport measurement^^. The corresponding inter¬ 
valley scattering length i v was obtained experimentally 
in Ref. [53] for CdsAs 2 . This length is connected to the 
scattering time via £ v = y/Dr v ^ where D = /x e k B T/e 
is the charge diffusion coefficient at temperature T. To¬ 
gether with mobility measurements, the inter-valley scat¬ 
tering time r v at T = 4K can be determined to be 
r v ~ 10 -9 s, similar to the theoretical estimate for NasBi. 


2. ARPES in finite magnetic fields 

ARPES experiments in magnetic fields are challeng¬ 
ing. External magnetic fields affect electron trajectories 
(especially those with low energies) and compromise an¬ 
gle resolution. Typical ARPES equipment is protected 
from external magnetic fields using /i-metal shields made 
from a metallic alloy with high magnetic permeability /i. 
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Dirac Semimetals 

Weyl Semimetal 

CdsAs 2 Ref. 

NasBi Ref. 

TaAs Ref. 

Re (cm 2 V 1 s 1 ) 

10 4 - 10 5 [HlfTTlfTsl and [Ml 

~ 10 3 -10 4 |5^ 

5 • 10 5 El 

r c (s) [Eq. dslf] 

~ 10“ 14 -10“ 13 

~ 10- 13 

~ 10- 12 

t v / t c [theoretical] 

- 

O 

T-1 

500 EH 

Ti / t c [theoretical] 

- 

10 3 eh 

Ti not defined 

r v (s) [experiment] 

O 

1 

CO 

- 

- 

«(t) 

7.6 • 10 5 El 

9.3 10 5 El 

1.1-10® H3 

1.3-10® HU 

1.5-10® El 

3.74 • 10 5 El 

3.1 - 3.6 • 10 5 (wi) EH 
2.6 - 4.3 • 10 5 (W 2 ) EH 

(f) 

~ 10 5 El 

3.3 10 5 El 

2.89 • 10 5 El 

3.4 • 10 5 (wi) EH 

4.1 • 10 4 (W 2 ) EH 


a Ilya Belopolski, private communication. 

b Measurement of £ v , the inter-valley relaxation time was obtained as described in the text. 


TABLE I. Summary of experimentally measured Fermi velocities v and v z , mobilities /x e , and theoretical estimates for the 
relevant scattering times. For TaAs W1 (W2) nomenclature classifies the eight (sixteen) Weyl points falling on (away) the 
k z = 27r/c plane. The scattering rate r c is estimated for an isotropic Weyl or Dirac cone with Ef = HvkF ~ lOmeV. 


This material offers the magnetic field lines a path with a 
low magnetic resistance (or reluctance) that is inversely 
proportional to /x, preventing them going to the energet¬ 
ically costly exterior. In typical ARPES experiments, for 
instance the set-up of Ref. |49l high resolution requires 
the field inside the ARPES lens to be of the order of 
rT or less. In practice up to rsj 6 /iT can be handled^. 
These small magnetic fields typically induce a rigid shift 
of the electron trajectories that can in principle be cor¬ 
rected^. As we discuss in the next section, already these 
small fields can result in observable values of the chiral 
potential difference 5r ~ 4meV. 

To allow for measurements in larger magnetic fields, 
without altering electron trajectories, it would be advan¬ 
tageous to go one step beyond existing set-ups. We sug¬ 
gest that a possible experimental design can be based on 
early i deas used to perform ARPES on ferromagnetic ma¬ 
terials 5 ^ 57 . Stray fields around ferromagnetic materials, 
like nickel, can severely alter photo-electron trajectories 
and jeopardize ARPES measurement accuracy. An ele¬ 
gant experimental solution studied these materials in the 
so-called picture frame geometry or remnant stat d 50 * 5 ^ 
(see Fig. |5|. In this geometry, the magnetic field lines 
are confined within the material in a magnetic loop, min¬ 
imizing their effect on electron trajectories that can spoil 
ARPES measurements. To reach a uniform field within 
a Weyl or Dirac semimetal sample, while minimizing the 
external stray fields, one possible solution is to open up 
a small gap in the magnetic circuit as shown in Fig. |5j 
In this set-up, a very uniform field is obtained within the 
gap, with possible additional fringing fields depicted as 
curved field lines. The effect of the fringing field is to 
increase the effective area of the gap, thus reducing its 
reluctance and increasing the effective magnetic field felt 


by the sample 58 . The magnitude and spatial extent of 
the fringing field can therefore be optimized by increasing 
the gap’s length or cross section. In addition, if the sam¬ 
ple cross section is designed to be slightly larger than the 
cross section of the gap, the fringing fields could enter the 
sample as well, not affecting photo-electron trajectories 
as shown in Fig. [5] Finally we note that external current 
loops were commonly used in experiments to orient the 
magnetization 50 5 ^. Therefore it is plausible that these 
set-ups can admit modifications to incorporate electric 
fields as well. 


3. Estimates of the chiral chemical potential 
difference 


To observe the note-shaped pattern that signals the 
chiral anomaly it is essential that the induced chemical 
potential difference Sr is larger than the energy resolu¬ 
tion of the experiments, which is of the order of ~ meV— 
if magnetic fields are kept of the order of ^ /iT. To esti¬ 
mate the difference Sr = rl — Rr, we use the definition 
of the chiral chemical potentials given by equation (3) of 
the main text. For the particular case where /i eq = 0, 
this difference is given by 


5/1 = 2 


h v 2 v z e 2 r v E • 



(S2) 


where v and v z are the anisotropic cone Fermi veloci¬ 
ties defined in the Methods section of the main text and 
reproduced in the Table [I] as measured by ARPES for 
NasBi and CdsAs 2 - From Eq. (S2) we observe that it is 
the magnitude of the Fermi velocities, rather than r v or 
B, that can most effectively enhance the chiral anomaly. 





























FIG. 5. Proposed magnetic circuit for studying the chiral 
anomaly and ARPES in a magnetic field, based on ARPES 
experiments on ferromagnetic materials^. The sample is in¬ 
cluded in a magnetic circuit to minimize stray fields effects on 
photo-electron trajectories. The yellow arrows show the con¬ 
fined magnetization in the ferromagnetic material (e.g. nickel 
or iron) shaped into the picture-frame geometry. The Weyl 
or Dirac semimetal sample, a photon of energy huo and a 
photo-electron with energy and momentum {E, k} are de¬ 
picted schematically as a blue rectangle, a purple curved ar¬ 
row and a red straight arrow respectively. 


To calculate S/i we choose two different magnetic field 
values |B| = 6 //T and |B| = lmT. The first corre¬ 
sponds to a reasonable magnetic environment that can 
be achieved in ARPES as discussed in the previous sec¬ 
tion, while the second is a very conservative estimate of 
the magnetic fields that can be achieved with the frame 
geometry also introduced above. All our following esti¬ 
mates are calculated for |E| = 1 • 10 1 2 3 4 5 6 7 Vm -1 . For the 


relaxation rates, we choose the experimentally obtained 
time r v = 10 -9 s for Cd 3 As 2 and the theoretical esti- 
mate r v /r c = 10 4 for NasBi, both justified in Sec. |A l| 
These are also conservative estimates since r v is expected 
to be even larger for small chemical potentials 4 ^. Given 
these values, we find for CdsAs 2 and field strengths of 
|B| = lmT (6/iT) a chiral chemical potential difference 
of Sfi ~ 30meV (5.4 meV). For NasBi, we find that 
5/i ~ 15meV (2.7meV) for |B| = lmT (6/iT). 


Appendix B: Case example: NasBi 

The main features that can lead to a direct visualiza¬ 
tion of the chiral anomaly in Weyl and Dirac semimetals 
via photoemission spectroscopy are presented in the main 
text. In order to make definite contact with current ex¬ 
perimental state of the art, we discuss the features of 
the chiral anomaly for NasBi. To this extent, we com¬ 
bine the experimental constants in Table [I] with param¬ 
eters obtained from ab initio calculations for this mate¬ 
rial that define the realistic low energy model described 
in Ref. |24j The equilibrium ARPES spectra computed 
for such model, as detailed in the main text, is shown 
in Fig. [ 6 ^a). The figure shows a clear pair of curved 
Fermi arcs that connect two Dirac nodes. Distinguishing 
clearly these pair of arcs is well within experimental res¬ 
olution, as confirmed by recent experiments [cf. Fig. 2A 
of Ref. mi- Upon applying external fields that satisfy 
E • B 7 ^ 0, the chiral anomaly induces a chiral chemical 
potential difference of the order S/a ~ lOmeV, as esti¬ 
mated in the previous section. In Fig. [ 6 ^b), we show 
the chiral anomaly induced ARPES spectrum for NasBi 
corresponding to such value of S/i. The Fermi arcs are 
partially occupied and thus end before reaching the bulk 
Fermi surface, the latter not visible within this resolution. 
The observation of such features would directly visualize 
the chiral anomaly and confirm its effect on the surface 
electronic spectrum predicted in the main text. 
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